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ABSTRACT 

We introduce the NEXUS algorithm for the identification of Cosmic Web environ- 
ments: clusters, filaments, walls and voids. This is a multiscale and automatic mor- 
phological analysis tool that identifies all the cosmic structures in a scale free way, 
without preference for a certain size or shape. We develop the NEXUS method to 
incorporate the density, tidal field, velocity divergence and velocity shear as tracers of 
the Cosmic Web. We also present the NEXUS+ procedure which, taking advantage 
of a novel filtering of the density in logarithmic space, is very successful at identifying 
the filament and wall environments in a robust and natural way. 

To asses the algorithms we apply them to an N-body simulation. We find that all 
methods correctly identify the most prominent filaments and walls, while there are 
differences in the detection of the more tenuous structures. In general, the structures 
traced by the density and tidal fields are clumpier and more rugged than those present 
in the velocity divergence and velocity shear fields. We find that the NEXUS+ method 
captures much better the filamentary and wall networks and is successful in detecting 
even the fainter structures. We also confirm the efficiency of our methods by examining 
the dark matter particle and halo distributions. 

Key words: Cosmology: theory - large-scale structure of Universe - Methods: data 
analysis - techniques: image processing 



1 INTRODUCTION 

Early attempts to map the large scale distribution of galax- 
ies in the universe (Gregory & Thompson 1978; Geller & 
Huchra 1989; de Lapparent et al. 1986; Shectman et al. 1996) 
revealed that galaxies are far from being evenly distributed 
across the nearby Universe. On the contrary, the mass dis- 
tribution delineated by galaxies seems to form an intricate 
network of compact and dense associations interconnected 
by tenuous "bridges" or "filaments" surrounded by surpris- 
ingly vast empty regions (Kirshner et al. 1981). Preliminary 
studies suggested that the universe on large scales could be 
described as a cellular system (Joeveer & Einasto 1978) or 
a Cosmic Web (Bond et al. 1996). This has been confirmed 
in recent times by large galaxy surveys such as the 2dFGRS 
(Colless & et. al. 2003), the Sloan Digital Sky Survey (e.g. 
Tegmark & et. al. 2004) and the 2MASS redshift survey 
(Huchra & et al. 2005). 

The Cosmic Web can be seen as the most prominent 
manifestation of the anisotropic nature of gravitational col- 
lapse, the motor behind the formation of structure in the 
cosmos (Peebles 1980). N-body computer simulations have 
illustrated how a primordial field of tiny Gaussian density 
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perturbations transforms into a pronounced and intricate 
filigree of filamentary features, dented by dense compact 
clumps at the nodes of the network (Jenkins & et al. 1998; 
Colberg et al. 2005; Springel & et. al. 2005; Dolag et al. 
2006). The description of the Megaparsec matter distribu- 
tion as an interconnected network or a cosmic web is not a 
coincidence. Even early computer simulations indicated the 
close connection between each morphological component, 
namely that clusters sit at the intersection of filaments and 
filaments are formed at the intersection of walls (Doroshke- 
vich et al. 1980; Melott 1983; Pauls & Melott 1995; Shapiro 
et al. 1983; Sathyaprakash et al. 1996). 

One of the main reasons for our interest in outlining 
the Cosmic Web concerns the question whether and how 
far the weblike environment influences the properties and 
evolution of galaxies. Recent iV-body simulations have found 
that the filamentary or sheetlike nature of the environment 
has a distinct influence on the shape and spin orientation of 
dark matter haloes (Aragon-Calvo et al. 2007a; Hahn et al. 
2007a,b; Paz et al. 2008; Hahn et al. 2009; Hahn 2009; Zhang 
et al. 2009). Other recent works (Jones et al. 2010; Tempel 
et al. 2012) have shown that indeed there is an alignment, 
even though weak, of galaxies and the filaments they lie 
within. In this paper we propose new robust and flexible 
methods that allow a better identification of the Cosmic 
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Web environments and hence help us to better understand 
how environments influence the formation and evolution of 
dark matter haloes and galaxies. 

1.1 Cosmic Web Detection 

Identifying the components of the Cosmic Web is a major 
challenge due to the overwhelming complexity of the indi- 
vidual structures as well as their connectivity, the lack of 
structural symmetries, its intrinsic multiscale nature and the 
wide range of densities found in the cosmic matter distribu- 
tion. Over the years, a variety of heuristic measures were 
forwarded to analyse specific aspects of the spatial patterns 
in the large scale Universe, but only recently these have lead 
to a more solid and well-defined machinery for identifying 
the Cosmic Web. Nearly without exception, these methods 
borrow extensively from other branches of science such as 
image processing, mathematical morphology, computational 
geometry and medical imaging. 

The connectedness of elongated supercluster structures 
in the cosmic matter distribution was first probed by 
means of percolation analysis, introduced and emphasized 
by Zel'dovich and coworkers (Zeldovich et al. 1982; Shan- 
darin & Zeldovich 1989; Shandarin et al. 2004; Shandarin 
et al. 2010), while a related graph-theoretical construct, the 
minimum spanning tree of the galaxy distribution, was ex- 
tensively analysed by Bhavsar and collaborators (Barrow 
et al. 1985; Graham & Clowes 1995; Colberg 2007) in an 
attempt to develop an objective measure of filamentarity. 
Both Colberg et al. (2005) and by Pimbblet (2005) set out 
to identify filaments and their adjoining clusters, using quite 
different techniques. 

More general filament finders have been put forward 
by a number of authors. Following specific physical criteria, 
Gonzalez & Padilla (2010) recently forwarded an interesting 
and promising combination of a tessellation-based density 
estimator and a dynamical binding energy criterion. A thor- 
ough mathematical nonparametric formalism involving the 
medial axis of a point cloud, as yet for 2-D point distribu- 
tions, was proposed by Genovese et al. (2010). It is based on 
a geometric representation of filaments as the medial axis of 
the data distribution. Also solidly rooted within a geomet- 
ric and mathematical context is the more generic geomet- 
ric inference formalism developed by Chazal et al. (2009). 
It allows the recovery of geometric and topological features 
of the supposedly underlying density field from a sampled 
point cloud on the basis of distance functions. Stoica et al. 
(2005, 2007, 2010) use a generalization of the classical Candy 
model to locate and catalogue filaments in galaxy surveys. 
This approach has the advantage that it works directly with 
the original point process and does not require the creation 
of a continuous density field. However, computationally it is 
very demanding. 

The more recent formalisms that are intent on charac- 
terizing the full range of weblike formalisms usually exploit 
the morphological information in the gradient and Hessian of 
the density field or potential field, i.e. the tidal field (see e.g. 
Aragon-Calvo et al. 2007a, b; Hahn et al. 2007a,b; Sousbie 
et al. 2008; Forero-Romero et al. 2009; Bond et al. 2010a,b). 
Morse theory (see Colombi et al. 2000) forms the basis of the 
skeleton analysis by Novikov et al. (2006) (2-D) and Sousbie 
et al. (2008) (3-D). It identifies morphological features with 



the maxima and saddle points in the density field and results 
in an elegant and mathematically rigorous tool for filament 
identification. However, it is computationally intensive, fo- 
cuses mostly on filaments and is strongly dependent on the 
smoothing scale of the density field. A more elaborate classi- 
fication scheme on the basis of the manifolds in the tidal field 
- involving all morphological features in the cosmic matter 
distribution - has been forwarded by Hahn et al. (2007a) 
(also see Hahn et al. 2007b; Forero-Romero et al. 2009). 

Instead of using the tidal field configuration, one may 
also try to link directly to the morphology of the density field 
itself. Though this allows a more detailed view of the intrica- 
cies of the multiscale matter distribution, it is usually more 
sensitive to noise and less directly coupled to the underlying 
dynamics of structure formation than the tidal field mor- 
phology. A single scale dissection of the density field into its 
various morphological components has been done by Bond 
et al. (2010a), and applied to N-body simulations and galaxy 
redshift samples (also see Bond et al. 2010b; Choi et al. 
2010). A more elaborate formalism is the Multiscale Mor- 
phology Filter (MMF), introduced by Aragon-Calvo et al. 
(2007b). It looks at structure from a scale-space point of 
view, treating the spatial structure in D dimensions in the 
context of an explicit D + 1 dimensional space (Florack et al. 
1992; Lindeberg 1998). The D + 1 dimensional scale space 
consists of the D-dimensional spatial structure at a range of 
spatial resolution scales. The MMF formalism subsequently 
assigns a local morphology based on an evaluation of the 
multiscale second order variations in the local density field. 
Instead of restricting the analysis to one particular scale, 
by evaluating the density field Hessian over a range of spa- 
tial scales and determining at which scales and locations 
the various morphological signatures are most prominent, 
the MMF explicitly addresses the multiscale nature of the 
cosmic structures. A somewhat similar multiscale approach 
was followed by the Metric Space Technique described by 
Wu et al. (2009) , who applied it to a morphological analysis 
of the DR5 of the SDSS. 

A more recent development is that of the Spineweb pro- 
cedure (Aragon-Calvo et al. 2010), which traces the various 
features of the cosmic web on pure topological grounds by in- 
voking the Watershed Transform (WT). The WT was intro- 
duced by Platen et al. (2007) as the basis of the Watershed 
Void Finder technique (Platen et al. 2007) which identifies 
cosmic voids with the watershed basins. The Spineweb pro- 
cedure elaborates on this, by identifying the central axis of 
filaments and the inner plane of walls with the boundaries 
between the watershed segments of the density field. While 
the basic Spineweb procedure involved one scale, the full 
procedure allows a multiscale topological characterization 
of the Cosmic Web Aragon-Calvo et al. (2010). However, 
to do so it must invoke some implicit assumptions on the 
connectivity of the various topological features. 

1.2 Intention and Outline 

The goal of this paper is to present two new algorithms 
(NEXUS and NEXUS+) for the detection of Cosmic Web 
environments and to assess their effectiveness. Elaborating 
on the multiscale scale-space context of the rudimentary 
density field Multiscale Morphology Filter (MMF, Aragon- 
Calvo (2007)), the NEXUS and NEXUS+ formalism repre- 
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sents a complete and versatile instrument for the structural 
and physical study of the Cosmic Web. Simultaneously tak- 
ing into account the multiscale nature of the cosmic mass 
distribution, NEXUS and NEXUS+ explicitly operate on a 
diversity of physical fields that are relevant to the formation 
and evolution of the Cosmic Web. The extension beyond 
the density field, towards the use of information contained 
in the tidal field, velocity divergence and velocity shear to 
trace the large scale structure, is a key aspect of NEXUS and 
NEXUS+. The new formalism allows us to compare the en- 
vironments traced by both the positional and velocity part 
of the phase space. We focus most of the analysis on the 
detection of filaments and walls, since these are the most 
challenging environments to identify. We find that our new 
methods are very efficient at tracing the Cosmic Web, re- 
sulting in very high quality filaments and walls. 

This study is organized as follows. In sections 2 and 3 
we describe the NEXUS and NEXUS+ methods which we 
use for the identification of the Cosmic Web, including a 
comparison of the two algorithms on a toy model. This is 
followed in section 4 with an extension of the tools to use 
a multitude of cosmological fields (density and tidal field 
versus the velocity divergence and velocity shear) as tracers 
of the cosmic environments. This way we take full advantage 
of the full 6-D information contained in phase space. 

The second part of the paper is focused on assessing 
how these methods cope with the complex and hierarchical 
structures present in the universe. To do so we use the DTFE 
density and velocity divergence from N-body simulations as 
inputs to our algorithms - see section 5. Section 6 presents 
the cluster, filament and wall environments identified in the 
simulation and compares the results of the different meth- 
ods. To asses the quality of the detections we look at the 
dark matter particle and halo distributions in each environ- 
ment and also study the effects of a multiscale versus single 
scale approach in tracing the Cosmic Web. Finally, we sum- 
marize our findings in section 8. 



2 NEXUS : GENERAL FORMALISM FOR 
MULTISCALE MORPHOLOGICAL 
ANALYSIS 

The NEXUS algorithm is a Scale-Space method for mor- 
phologically segmenting the Cosmic Web into its three dis- 
tinct features: clusters, filaments and walls. The environ- 
ment identification is performed in a scale and user indepen- 
dent way to account for the multiscale nature of the large 
scale structure, which is the result of the hierarchical evolu- 
tion of the cosmic mass distribution. The method is derived 
from the field of medical imaging (Frangi et al. 1998; Sato 
et al. 1998; Li et al. 2003) where it is used to identify nod- 
ules, vessels and walls in two- and three-dimensional images. 
An earlier and simpler version of the method was introduced 
in cosmology under the name Multiscale Morphology Filter 
(MMF) in Aragon-Calvo et al. (2007b). The MMF involved a 
rudimentary treatment of scale-space analysis and restricted 
itself to the use of the density field as tracer of the Cosmic 
web environment. 

In general, the scale-space formalism can be applied to 
any input field to detect point-, line- and sheetlike struc- 
tures in the field values, and thus lends itself to applications 



involving a range of quantities dynamically relevant for the 
formation and evolution of the Cosmic Web. Following this 
observation, we have embedded the scale-space formalism in 
the physical framework of Cosmic Web formation. 

For this purpose we have defined two classes of the algo- 
rithm, the NEXUS and the NEXUS+ formalism. The main 
difference between NEXUS and NEXUS+ concerns the fil- 
ter used for constructing the representation of the field at 
different resolutions in scale-space. The NEXUS technique 
uses a Gaussian filter for smoothing while NEXUS+ uses 
a Log-Gaussian filter (more on that in section 3). We will 
demonstrate in this study that they yield a substantially 
more realistic and robust representation of filaments, walls 
and their mutual connectivity, over the range of scales cov- 
ered by the scale-space representation. 



2.1 NEXUS : general algorithm description 

The NEXUS algorithm detects the point-, line- and sheet- 
like structures for a generic input field /. For large scale 
structure, these features correspond to clusters, filaments 
and walls. To keep the notations clear, we limit our discus- 
sion to the Cosmic Web environments, but there is no loss of 
generality. The NEXUS algorithm consists of the following 
six steps: 

(I) Applying a Gaussian filter of width R n to the input 
field. 

(II) Computing the Hessian matrix eigenvalues for the fil- 
tered field. 

(III) Assigning to each point a cluster, filament and wall 
signature using the Hessian eigenvalues . 

(IV) Repeating steps (I) to (III) over a range of smoothing 
scales (Ro, Ri, ■ ■, Rn), to construct the scale-space 
representation of the field. 

(V) Combining the results of all scales to obtain a scale 

independent cluster, filament and wall signature. 
(VI) Using physical criteria to determine the detection 
threshold corresponding to valid environments. 

In the following we elaborate on each step of the algorithm 
and give the details necessary for the implementation of the 
method. 



2.1.1 Step I: Applying Gaussian smoothing 

A Gaussian filter of width R n is applied to the input field 
/(x). This gives rise to a smoothed field /r„(x) given by: 

. . f d A k -k 2 R 2 II ?/, n ikx /,\ 

where /(k) is the Fourier transform of the input field /(x). 



1 The algorithm given here applies to detecting the point-, line- 
and sheet-like structures corresponding to maxima in the field val- 
ues. If we are interested in the same structures but for the minima 
of the field values, than we need to apply the same algorithm to 
— /. For example the Cosmic Web environments correspond to 
maxima in density but to minima in velocity divergence. 
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Figure 1. The amplitude of the kernel function Hij for a smooth- 
ing radius of 1 h~ 1 Mpc. The values were normalized such that 
Hij has a maximum value of 1. Note the logarithmic axes. 

2.1.2 Step II: Computing Hessian eigenvalues 
The Hessian of the filtered field is computed as: 

, 2 a 2 /fl„(x) 



Hij,R n (x) = R n 



dxidxj 



(2) 



where Hij t R n represents the i,j entry of the Hu n Hessian 
matrix. The R n term is a renormalization factor that has to 
do with the multiscale nature of the NEXUS algorithm. It 
makes sure that the same weight is assigned when compar- 
ing the Hessian value at different scales. Using Eq. (1), the 
Fourier transform of the Hessian reads: 



with H the Hessian kernel function given by 



(3) 



(4) 



The kernel function characterizes which Fourier components 
of the input field give contributions to the Hessian matrix. 
The dependence of the Hessian kernel on k is shown in Fig- 
ure 1. At a given smoothing scale R n , the H kernel has a 
peak at fc pca k = \/2/ R n , with a sharp drop-off for higher k 
and a linear fall for smaller k. Therefore, for a given scale 
R n , only the Fourier components of / around the peak fc pe ak 
will give an important contribution to the Hessian matrix. 

The NEXUS formalism depends only on the eigenvalues 
of the Hessian matrix, eigenvalues given by: 

det(H R Jx) ~ A„,ij n (x)) = 0, with Ai s= A 2 < A 3 . (5) 



2.1.3 Step III: Computing environment signature 

The eigenvalues of the Hessian matrix can be used to assign 
a cluster, filament and wall characteristic to every point x 
using the expected behaviour given in Table 1. This is the en- 
vironment signature and is denoted with 5(x). The first step 
in computing the signature is to define the shape strength 
T. This gives a quantitative description of the approximate 
relations given in the middle column of Table 1. The shape 



Table 1. Hessian eigenvalue relationships for the different envi- 
ronments of the Cosmic Web. The second column gives the qual- 
itative relationships between the eigenvalues (conditions that are 
implemented analytically in Eq. (6) ) while the third column gives 
strict eigenvalues constraints implemented in Eq. (7). 



Structure 


Soft constraints 


Strict constraints 


cluster 


[Ai| ~ |A 2 | ~ |A 3 | 


Ai < 0; A 2 < 0; A 3 < 


filament 


|Ai| ~ |A 2 | » |A 3 | 


Aj < 0; A 2 < 


wall 


|Ai| » |A 2 |; |Ai| » |A 3 | 


Ai < 



strength is defined as: 



A3 

I A l 

I As; 

e(i 



0(i- 



A 9 

Ai 



e i 



Us 
Al 



cluster 

filament 

wall 



(6) 



where we use the notation O(x) = x9(x) for clarity, with 
9(x) the step function (9(x) = 1 if x ^ 0, otherwise). 
The strength X is large when the eigenvalues at x corre- 
spond to a prominent structure and small otherwise. The 
cluster/filament/wall signature is defined as: 




(-Ai)0(-A 2 )0(-A 3 ) 

(-Ai)0(-A 2 ) 

(-Ai) 



cluster 

filament (7) 
wall, 



where the 0(—\ a ) factors (with a — 1,2,3) incorporate the 
right most column of Table 1. The |A a | term gives the in- 
tensity of the morphological feature and can be used to dis- 
criminate between real signals (large |A a |) and noise (small 
|A.|). 

2.1.4 Step IV: Computing the environmental signature 
over a range of smoothing scales 

The previous three steps are repeated over a range of 
smoothing scales (Ro, Ri, ■■, Rn). The hierarchy of smooth- 
ing scales is taken as R n = (\/2) n i?o with Ro the smallest 
scale at which one expects to find structures (Sato et al. 
1998). Taking an even smaller step between any two suc- 
cessive smoothing scales makes only minor differences. In 
practice we choose Ro equal to the grid spacing of the input 
field. We found that for the detection of the most promi- 
nent features of the Cosmic Web it is sufficient to consider 
smoothing scales in the range 0.5 /i -1 Mpc to 4 /i~ 1 Mpc. 
In this respect it is good to note that outstanding features 
of the Cosmic Web are visible within a particular range of 
scales centred around the transition scale between linear and 
non-linear structures. 

The result of this step is a signature function for 
each scale 5_r„ (x) which characterizes the environmental re- 
sponse of point x at the R n smoothing scale. 



2.1.5 Step V: Scale-space stacking 

The signature of the given set of scales is combined to ob- 
tain the overall signature. This is a scale independent map 
characterizing the degree to which the point x is part of a 
cluster, filament or wall. A structure of a given size will give 
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Figure 2. A 1 h~ 1 Mpc slice through: a) density field, b) node 
signature, c) filament signature and d) wall signature obtained 
from an N-body simulation. The white, orange and black show 
the high, medium and low values of density and environment 
signature. The green line contours show the regions with high 
node signature (panel c) and with high filament signature (panel 
d). 



the largest signature for a smoothing scale of the same size. 
Therefore, the overall signature at a point is the maximum 
signature over all the scales: 



<S(x) = max Sr„(x). (8) 

levels n 



2.1.6 Step VI: Computing the detection threshold 

The signature has a wide range of values, with the large 
one corresponding to strong structures and the small ones 
coming from noise and null detections. This can be appreci- 
ated in Figure 2, which shows the cluster, filament and wall 
signature. Therefore, the last step in the algorithm involves 
the use of physical criteria to find the threshold signature 




6 7 
cluster signature S r 




0.1 1 10 

filament/wall signature S^ w 

Figure 3. Upper panel: The dependence of the fraction of clus- 
ters with density larger than the virial density versus the cluster 
signature iS c . The intersection of the gray lines shows the cluster 
detection threshold. Lower panel: Determination of the detection 
threshold for filaments and walls. The peak of AM 2 (shown by 
the gray vertical lines) corresponds to the signature threshold for 
filament and wall identification (see text for details). 



that discriminates between valid and invalid detections. Sig- 
nature values larger than the threshold correspond to real 
structures while the rest are spurious detections. The thresh- 
old signature for clusters is found by requiring that the iden- 
tified objects are virialized, whereas for filaments and walls 
the threshold is given by the dependence of the filament / wall 
mass with environmental signature. 

The procedure to determine the signature threshold for 
cluster detection is illustrated in the upper panel of Figure 
3. Clusters are the largest and most recently formed fully 
virialized objects (Voit 2005). We use this definition to de- 
termine the signature threshold for cluster identification. We 
test for virialization by requiring that the average density of 
the cluster is larger than A = 370, which is the value given 
by the spherical collapse model at z = (Gunn & Gott 
1972). From Figure 3 it can be seen that the fraction of ob- 
jects with an average density larger than the virialization 
threshold changes very fast from to 1 as we increase the 
cluster signature S c . We then take the signature threshold 
as the value where half of the objects have a density larger 
than A. 

The filament and wall identification is performed by lim- 
iting our detections to only the most prominent filamentary 
and wall regions. We find that the same method can be sue- 
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cessfully used for the recognition of both filaments and walls. 
Let us denote with Mf(Sf) the mass in filaments with a sig- 
nature value larger or equal to <S/. As Sf decreases, more 
and more regions are included and hence Mf(Sf) increases. 
Most of the change in this function is restricted to a small 
range in Sf values and it gives a natural way of discriminat- 
ing between real and spurious detections. In Appendix A we 
show that the mass change with signature: 



a) Input image 



AM 



dM 2 



d log Sf 



(9) 



gives a natural and robust method of defining the most 
prominent filamentary components of the Cosmic Web. Sim- 
ilarly, in the case of walls, we can define AMI 2 , using the 
above equation with Mf and Sf replaced by their corre- 
sponding quantities for wall environments, M w and S w . The 
quantity M W (S W ) is the mass in regions that have a wall sig- 
nature value larger or equal to S m . The AM 2 dependence 
for both filaments and walls is shown in Figure 3. We use the 
pronounced AM 2 peak to delineate the valid environments, 
which are the points with signatures larger than the position 
of the AM 2 peak. All other points with smaller signature 
are considered null detections. This threshold method repro- 
duces very well the filamentary and wall network visible in 
both the cosmic density and velocity divergence fields. 

The algorithm performs the environment detection by 
applying the above steps first to clusters, then to filaments 
and finally to walls. This sequence (first clusters, then fila- 
ments and finally walls) has to be followed due to presence 
of anisotropic clusters and filaments that give mixed envi- 
ronmental signatures. This can be appreciated from panel c 
of Figure 2 where on top of the filament signature we show 
the contours corresponding to large cluster signature. We 
see that there are several regions that have both a large 
cluster and filamentary characteristic. This is due to non- 
spherical clusters that have a large filamentary signature. 
Similar anisotropic cluster/filaments may give a strong wall 
signature (see panel d of Figure 2). To overcome this cross- 
contamination, a point is part of a filament only if it was 
not previously identified as in a cluster. Similarly a point is 
in a wall if it was not previously identified as part of a clus- 
ter or filament. This procedure makes sure that each point 
is assigned a single classification: cluster, filament, wall or 
field (everything else that is not a cluster, filament or wall). 



3 NEXUS+ : LOGARITHMIC FORMALISM 
FOR MULTISCALE MORPHOLOGICAL 
ANALYSIS 

The NEXUS algorithm is very efficient in detecting the en- 
vironments in a field / where all the structures correspond 
to the same order of magnitude values of /. However, the 
method faces some challenges when the structures in / are 
present over orders of magnitude in field values. To better 
understand this, we present a test example in Figure 4. It 
shows three filaments characterized by different intensities: 
1, 10 and 100 (from left to right). For the NEXUS method 
to identify all the three filaments, the threshold needs to be 
so low that the stronger filaments are detected as extending 
much beyond their input data boundaries. The higher inten- 
sity peaks give a significant signal even at large distances, 




Figure 4. The filaments detected using the NEXUS and 
NEXUS+ methods when applied to a test image. The input test 
image (upper panel) contains three filaments with the same width 
but of different intensities: 1, 10 and 100 (from left to right). 
Frame b) shows the NEXUS filament signature with a thresh- 
old low enough such that also the weakest filament (left most 
one) is visible. Panel c) depicts the same as frame b) but for the 
NEXUS+ method. 



due to the combination of the Gaussian filter not dropping 
off fast enough and the high field value of the peak. 

One way to remedy this problem is to replace the Gaus- 
sian filter with a new smoothing method that takes into ac- 
count the large range in values. For that we introduce the 
Log-Gaussian filter, which is a Gaussian filter in logarithm 
space. By replacing the Gaussian filter in NEXUS with the 
Log-Gaussian filter we obtain the NEXUS+ formalism. The 
results of the new method are presented in frame c of Figure 
4. It clearly shows that the new filtering procedure recovers 
much better the three filaments. 



3.1 NEXUS+ : general algorithm description 

The main difference between NEXUS and NEXUS+ is the 
use of the Log-Gaussian filter instead of the Gaussian one. 
The steps of the NEXUS+ algorithm are the same as the 
steps of NEXUS with the exception of steps (I) and (II). 



3.1.1 Step I: Applying Log-Gaussian smoothing 

A Log-Gaussian filter of width R n is applied to the input 
field /. To this end, we introduce the field g, the logarithm 
of field /, 

9 = log 10 / , (10) 
and the field gn n , the smoothed logarithm at scale R„, 



9R„ 



dyg(y) W ,« B (x,y), 



(11) 



with Wg,r„ the Gaussian filter of width R n . 

Following the introduction of these quantities, the ap- 
plication of the Log- Gaussian filter consists of three main 
steps, 
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1. Computing the logarithm of the field f,g = log 10 /. 

2. Applying the Gaussian filter of width R n to g to ob- 
tain the smoothed logarithm gn n . 

3. Computing the smoothed field fn n by taking the ex- 
ponential of the smoothed logarithm gn n . 

In practice, we perform the convolution of the field g with 
the Gaussian filter in Fourier space, 

d 3 k 



9R n (x) = 



(2tt) 3 



5(k) f 



(12) 



involving the simple multiplication of the Fourier field com- 
ponent g(k) with the Gaussian exponential, 



9R n ( k ) = e 



-k 2 Rfj2 . 



S(k)- 



(13) 



Subsequently, the resulting NEXUS+ smoothed field fa n is 
obtained by evaluating, 



/fl„(x) = Cjj„ 10 9 



(14) 



The variable Cr^ is a multiplication constant that assures 
the mean of the input field is the same before and after 
filtering. 

3.1.2 Step II: Computing Hessian eigenvalues 

The second step is the same as for the NEXUS algorithm, 
but since the smoothing filter is different some of the equa- 
tions will also change. Now the Hessian of the smoothed 
field: 

, 2 a 2 /«„(x) 



Hij,R n (x) = R n 



dxida 



can be written in Fourier space using: 

Hij,R n (k) = -kikjRlfRjk). 



(15) 



(16) 



Please note that now one cannot formulate /fj n (k) as a sim- 
ple analytical expression as in the case of the NEXUS algo- 
rithm's Eq. (3). For the NEXUS+ algorithm one needs to 
perform steps (I) and (II) separately. 

The rest of the steps are the same as NEXUS steps (III) 
to (VI) described in section 2.1. It is important to note that 
because NEXUS+ uses the logarithm of / it can only be 
applied to input fields that have positive values at every 
point. 

3.2 NEXUS+ on the density field 

The major challenge of structure detection lies in the fact 
that the nonlinear density field, following its evolution, 
ranges over many orders of magnitude between the under- 
dense and overdense regions. Structures are present over the 
whole range of values in density. To deal with this challenge 
we can use two approaches: either take the density logarithm 
(see section 4.3) or use a different algorithm that takes into 
account the approximative log-normal shape of the density 
distribution. Here we take the former approach and apply 
the NEXUS+ algorithm to the density field. 

The strength of the NEXUS+ algorithm can be easily 
appreciated if one compares the density field using Gaussian 
versus Log-Gaussian smoothing. While there is a one to one 
mapping between density and density logarithm, this rela- 
tion does not hold when one compares the smoothed density 
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Figure 5. Comparing Gaussian and Log-Gaussian density 
smoothing (see text for details). The upper panel shows the initial 
density field, while the central and lower panels give the Gaussian 
and respectively Log-Gaussian smoothed density. Both cases were 
obtained using a 1 1 Mpc smoothing. The scale in the lower two 
panels was selected by fitting the density histogram with a log- 
normal distribution and plotting the values in the peak — 3cr to 
peak+3cr range (with peak and <r the peak and standard deviation 
of the log-normal distribution). 
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with the smoothed logarithm of the density. The results of 
the two methods are shown in Figure 5. We immediately 
observe that the Gaussian filtered density is dominated by 
several peaks with a typical spherical shape. On the other 
hand, the Log-Gaussian results seem to trace much better 
the large scale structure. 

Most of the differences between the two results come 
from the very high density peaks. When applying a Gaus- 
sian smoothing, these higher density peaks gets smoothed 
up to large distances and dominate the signal coming from 
other less dense neighbouring regions. This leads to a loss 
of information about the large scale structure around these 
peaks. In the case of Log-Gaussian smoothing, by taking the 
logarithm of the input field the contrast between these very 
high density peaks and their neighbourhoods is greatly re- 
duced. Therefore the contribution of the peaks will not be 
dominant, even though the Log-Gaussian filter has the same 
spatial extension as the Gaussian one. 



4 TRACER FIELDS OF THE COSMIC WEB: 
EXTENDING THE NEXUS ALGORITHM 
BEYOND DENSITY 

There are various methods that attempt to identify the com- 
ponents of the Cosmic Web. These not only implement dif- 
ferent detection techniques, but in many cases differ in the 
nature of the field used to trace the underlying cosmic struc- 
ture. In other words, the variation in the results of different 
methods should not only be ascribed to the algorithms used, 
but also to the differences in the tracer fields. In this sec- 
tion we extend the NEXUS method to a multitude of Cos- 
mic Web tracers: density, tidal field, velocity divergence and 
shear as well as to the density logarithm. By doing so we 
not only find the field with the best footprint of the Cosmic 
Web, but also gain better understanding of the evolution 
and structure of the Cosmic Web. 

The most widely used tracers of the Cosmic Web is 
the density (Aragon-Calvo et al. 2007b; Sousbie et al. 2008, 
201 1) and the tidal field (Hahn et al. 2007a; Forero- Romero 
et al. 2009). On the other hand, the use of the density log- 
arithm for the classification of the Cosmic Web is a novel 
method that we introduce here. Later on we will argue why 
this is a natural structure tracer field that one should con- 
sider. While there have been a fair share of methods using 
the positional information of the phase space, there are very 
few works that use the velocity field as Cosmic Web tracer. 
Most interesting is the work by Shandarin (201 1), which 
emphasized the importance of the velocity field for under- 
standing the emerging patterns in the matter distribution. 
Based on this, he used the variance of the velocity field as a 
measure of the local environment. In a sequence of studies 
that follows up on this idea (Abel et al. 2011; Shandarin 
et al. 2012; Neyrinck 2012), the full phase space structure 
of the mass distribution is used for an impressively accurate 
dynamical characterization of morphological structure of the 
Cosmic Web. Following this promising avenue, we are also 
working on relating the NEXUS and NEXUS+ formalism to 
the structure found by these phase space based methods. 



Table 2. The methods resulting from the extension of the 
NEXUS algorithm to several tracer fields of the Cosmic Web. 



Method name 


Tracer field 


NEXUS.den 


density field 8 


NEXUS_tidal 


tidal field T 


NEXUS.denlog 


density logarithm field log 10 (l + S) 


NEXUS_veldiv 


velocity divergence field 6 


NEXUS.velshear 


velocity shear field a 



4.1 NEXUS.den: tracing the Cosmic Web using 
the density field 

The density is one of the obvious fields used for environmen- 
tal detection due to the sharp contrast of the clusters and 
filaments compared to most of remaining void dominated 
volume. To apply the NEXUS algorithm on the density one 
has to just insert the density S in Eq. (3). For simplicity, we 
denoted this method as NEXUS.den . To better understand 
the behaviour of the NEXUS_den method we need to rewrite 
the Hessian matrix (see Eqs. (3)-(4)) to: 

= ~^(k) k 2 Rl e- k *<l\ (17) 

where the initial equation was multiplied by the unit factor 
^-fe 2 . It is immediately obvious that the Hessian is given by 
two distinct parts: the tidal field 2 multiplied by a band-pass 
filter. The band pass filter is made of two distinct compo- 
nents: the fc 2 -R 2 high pass filter and the Gaussian e~ h R ™/ 2 
low pass filter. Simple calculations show that the maximum 
of the band pass filter is at k = while the shape of the fil- 
ter is very similar to the one in Figure 1 . Therefore detecting 
the Cosmic Web structures in the density field is equivalent 
to identifying those structures in a band pass filtered tidal 
field. 

When applying the NEXUS-den formalism, an addi- 
tional step has to be taken on top of those described in 
section 2.1 and mask the density field when detecting fila- 
ments and walls. For filaments identification we need to set 
the density to in the cluster regions. In the absence of 
this mask, the cluster regions will give a large, unrealistic, 
filamentary signature 3 . Similarly when identifying walls, we 
need to set the density to in both the cluster and filament 
regions. 



2 The first part of Eq. (17) is the same as the Fourier transform 
of the tidal field given by Eq. (21) in section 4.2 up to the multi- 
plication factor AirGp. This factor has no effect on the final results 
since it only rescales the Hessian eigenvalues. 

3 The filamentary signature of cluster regions will be large, even 
though the filamentary shape strength Ifli am ent given by Eq. (6) is 
small in those regions. This is since the filamentary signature (see 
Eq. (7)) depends on | A2 1 which has very large values in the cluster 
regions and will compensate for the small values of Xfliament- So 
this additional | A2 1 factor that discriminates between signal and 
noise also introduces false detections. More generally, this false 
detection problem is important when the typical values of the 
tracer field for the different environments are orders of magnitude 
apart. It can be easily corrected by using the masking procedure 
described in the text. 
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Figure 6. The Fourier transform amplitude for the density, den- 
sity logarithm, gravitational potential, velocity divergence and ve- 
locity potential fields. The spectra were obtained directly from the 
DTFE interpolated fields used as input for the NEXUS method. 
The curves were shifted vertically to better emphasize the differ- 
ences between the two gray lines which mark the peak of the "Ha 
function with R = 4 (left line) and R = 0.5 h _1 Mpc (right line). 
The two smoothing radii represent the upper and lower limits of 
the smoothing scales set used in the NEXUS algorithm. 



4.2 NEXUS tidal: tracing the Cosmic Web using 
the tidal field 

The tidal field is the driver of anisotropic gravitational col- 
lapse and it is an essential ingredient for the formation and 
evolution of the Cosmic Web (Zel'dovich 1970; Gurbatov 
ct al. 1989; Bond et al. 1996). It is only natural to use it for 
the detection and understanding of the cosmic structures 
(Hahn et al. 2007a; Forero-Romero et al. 2009). This mode 
of the NEXUS method is indicated as NEXUS_tidal. 
The tidal field is given by 



Ty(x) 



d 2 0grav(x) 

dxidxj 



(18) 



with <f)grav the gravitational potential. The latter is related 
to the density via the Poisson equation: 



(x) = 4ttGp<5(x), 



(19) 



where G is the gravitational constant, p is the average matter 
density and l+8(x) = p(x) /p is the overdensity. The Poisson 
equation is easily solved in Fourier space to obtain: 



,(k) 



(20) 



leading to the following expression for the Fourier compo- 
nents of the tidal field: 

_ kj. k 7 



T i3 (k)=47rG^<5(k). 



(21) 



Identifying the cosmic environments traced by the tidal 
field is done by applying the NEXUS algorithm on the gravi- 
tational potential. The Hessian matrix of the potential <^ gra v 
given by Eq. (20) reduces to: 



k ; 2 -R„/2 



(22) 



This is exactly the tidal field smoothed over with a Gaussian 
filter. 




X (Mpc/h) 



Figure 7. The filament signature computed using: a) 
NEXUS_den and b) NEXUS_tidal methods. The graph shows a 
thin 1 h -1 Mpc slice. The white, orange and black correspond to 
high, medium and low signature values. 



The difference between the Hessian matrix of 
NEXUS_tidal and NEXUS_den consists in the additional 
k 2 R„ high pass filter present in the case of the second 
method - compare Eqs. (17) and (22). Therefore variations 
in the result of the two methods comes from excluding the 
low frequency modes in the NEXUS_den case and not in 
the NEXUS_tidal . The same conclusion can be reached 
by looking at the Fourier transform amplitude of the input 
fields: the density versus the gravitational potential. These 
are shown in Figure 6. The first is more flat, while for the sec- 
ond the low frequencies have much larger amplitudes. This 
means that the large scale modes give a much larger contri- 
bution for NEXUS.tidal than for NEXUS.den. These effects 
are illustrated in Figure 7. The NEXUS_den environments 
have a very clumpy appearance and are very sensitive to 
small scale structures. On the other hand, NEXUS_tidal is 
only responsive to the large scale modes and cannot trace 
the smaller details of the matter distribution. 

It is important to note that the environment character- 
istics in the gravitational potential are different from the 
ones in the density field. According to the Cosmic Web 
theory, the clusters, filaments and walls are given by the 
strength and sign of the first, second and third eigenvalues 
of the tidal field (Bond et al. 1996). This can be easily im- 
plemented within the NEXUS framework by changing the 
environmental signature from Eq. (7) to: 



(23) 



with a = 1 for clusters, a = 2 for filaments and a — 3 for 
walls. 

As in the case of the NEXUS_den method, we need 
to apply a cluster mask when identifying filaments and a 
combined cluster and filament mask when identifying walls. 
This procedure sets the density to in the mask regions, af- 
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ter which the gravitational potential is computed using Eq. 
(19). 

4.3 NEXUS_denlog: tracing the Cosmic Web 
using the density logarithm 

Using the density logarithm for Cosmic Web detection is an 
approach that has not been explored until now. We were 
motivated to apply the NEXUS algorithm to the density 
logarithm and not the density itself because of a multitude 
of reasons: 

• The NEXUS method works best when all structures 
correspond to similar values in the input field, while the den- 
sity ranges from 0.01 in underdense versus 10 4 and higher 
in overdense regions. We expect to find structures over or- 
ders of magnitude in density values and simply using density 
biases the results towards high density structures. By tak- 
ing the density logarithm the orders of magnitude difference 
is reduced to values of —2 in voids to around 4 in cluster 
regions 4 . 

• The non-linear density field is close to a lognormal dis- 
tribution, when smoothed on scales of a few Mpc (Coles & 
Jones 1991). 

• The large scale structure is best made visible when ren- 
dering the density logarithm and not the density itself. 

The NEXUS_denlog method consists in replacing the in- 
put field / in the NEXUS algorithm with log 10 (l + 5). The 
main difference between this method and NEXUS_den con- 
sists in the reduced contrast between underdense and over- 
dense regions as well as a much steeper spectrum to- 
wards large scales for the density logarithm (see Figure 
6). Because of the reduced contrast between underdense 
and overdense regions when looking at the density loga- 
rithm, there is no need to apply the mask described for the 
NEXUS_den method. 

While we were motivated by the same reasons as above 
to develop the NEXUS+ algorithm, there is a large dif- 
ference between the NEXUS+ and NEXUS_denlog meth- 
ods as will be clearly visible in the results of section 
6. In NEXUS_denlog we identify the Cosmic Web using 
the logarithm of the density log 10 (l + 8), while for the 
NEXUS+ method we trace the environments using the den- 
sity field smoothed with the Log-Gaussian filter. 

4.4 NEXUS_veldiv: tracing the Cosmic Web using 
the velocity divergence 

The NEXUS_den, NEXUS_tidal and NEXUS_denlog meth- 
ods use only half of the phase space, the positional infor- 
mation, for identifying the elements of the Cosmic Web. It 
is interesting from both a theoretical and practical point of 
view to see how the remaining phase space can also be used 
to trace the large scale structure. The natural candidates for 
this are the velocity divergence and the velocity shear, due 
to the one-to-one connection between these quantities and 
the density and tidal field in the linear regime. 

4 Note that the DTFE density field will always have a density 
different from zero even in the emptiest voids. The typical DTFE 
density contrast in voids at redshift z = is 0.01 to 0.1. 



The velocity divergence is defined as: 

<?( x ) = iv-v(x) (24) 

where we divide by the Hubble factor H such that 9 is a 
unitless quantity. The velocity divergence is easily computed 
as an output of the DTFE method (Bernardeau & van de 
Weygaert 1996; Romano-Diaz & van de Weygaert 2007). 
According to linear theory, the velocity divergence is related 
to the density field via: 

0(x) = -/5(x), (25) 

with / the linear velocity growth factor (see Peebles 1980). 
So in the linear regime, any structure in the density field 
should also be present in velocity divergence. A similar re- 
lation between 9 and 5 holds true also for the non-linear 
regime, but in a more complex way (for details see Nusser 
et al. 1991; Chodorowski & Lokas 1997; Bernardeau et al. 
1999). The differences between the structures detected us- 
ing density versus velocity divergence probe the effects of 
the non-linear evolution on the Cosmic Web components. 

The NEXUS_veldiv method is the application of the 
NEXUS algorithm on the negative of the velocity diver- 
gence —9. We choose the minus sign because of Eq. (25). 
The NEXUS_veldiv Hessian matrix is given by: 

H ijlRn (M) = ^m k 2 Rl e~ k *<l\ (26) 

This is the product of the velocity shear given by Eq. (32) 
multiplied by a band pass filter. This is exactly the same 
as for the NEXUS_den method, but with the tidal field re- 
placed by the velocity shear. The main difference between 
NEXUS_veldiv and NEXUS_den can be easily seen in Fig- 
ure 6: the velocity divergence and density have the same 
Fourier components at large scales, but the density has a 
more flattened drop at smaller scales. 

In contrast to the NEXUS_den method, we choose not 
to apply a mask for the velocity related methods. While 
there are still orders of magnitude variation in the velocity 
divergence between overdense and underdense regions, this 
difference is not as large as for the density field. The ma- 
jor challenge in applying a mask arises because the velocity 
divergence can take both positive and negative values, so 
there is no a priori well motivated value that we can use in 
the mask regions. 

4.5 NEXUS_velshear: tracing the Cosmic Web 
using the velocity shear 

The velocity shear is the symmetric part of the velocity gra- 
dient, with the ij component defined as: 

, , if dvj dvi \ ,„„. 

where Vi is the i component of the velocity. We normalize 
the velocity shear by the Hubble constant to keep the same 
notations as in the case of the velocity divergence. To obtain 
the velocity shear, the velocity is rewritten as the sum of the 
potential and rotational flows: 

v = V^voi + V x A vc i (28) 

where (f> vc \ is the scalar velocity potential and A vc i is the 
vector potential. Inserting this last equation into the velocity 
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shear gives: 



100 



,(x) = 



1 <9 2 0vcl(x) 

H dxidxj 



(29) 



Therefore, the velocity shear depends only on the velocity 
potential. It is interesting to notice that the velocity diver- 
gence is also given only by the velocity potential via: 



— V 2 VC 1. 

tl 



(30) 



This last equation can be inverted and used to solve for the 
potential </> ve i, to obtain its Fourier components as: 

1 



,(k) = -H-9(k). 



(31) 



On the basis of this equation, one can infer the velocity po- 
tential starting from the velocity divergence output of the 
DTFE method. Combining these, we obtain that the com- 
ponents of the velocity shear can be expressed in terms of 8 
as: 



(32) 



Thus the relation between the velocity shear and the tidal 
tensor in the linear regime is given by: 



(33) 



The NEXUS_velshear method involves using the neg- 
ative of the velocity potential — (f> vcl cLS input field to the 
NEXUS algorithm. The negative sign comes, as in the case 
of NEXUS_veldiv, from the minus in relation Eq. (33). In- 
serting the expression for <j) vc \ in the Hessian matrix gives: 



H 



ij,R. 



.(k) = -^0(k) 



(34) 



As in the case of the NEXUS_tidal method, the environment 
characteristics are different in velocity shear compared to 
velocity divergence. So the environment signature has to be 
changed to the expression given by Eq. (23). 
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Figure 8. A 1 h~ Mpc slice from the N-body simulation illus- 
trating the DTFE density 1 + <5 (upper panel) and absolute value 
of the velocity divergence 9 (lower panel). 



5 N-BODY SIMULATIONS AND HALO 
CATALOGUES 

To test our structure finding algorithms, we apply them to 
cosmological N-body simulations containing only dark mat- 
ter particles. We adopted the ACDM cosmological model 
with n m = 0.26, Q A = 0.74, h = 0.71, a 8 = 0.8 and 
n s — 1. We performed two 512 3 particle simulations in a 
100 h~ 1 Mpc and 200 ft - Mpc periodic boxes. The force res- 
olution was fixed in comoving coordinates up to z = 4, to 
15 h~ 1 kpc and 21 /i -1 kpc. Afterwards it was fixed in phys- 
ical coordinates to 5 h~ kpc and 7 h~ kpc, respectively. 

The 100 h~ Mpc simulation (mass resolution of 5.4 x 
10 8 /i _1 Mq) was chosen in order to resolve haloes to a few 
times 10 10 /i -1 Mq and at the same time to have a reasonable 
cosmological volume whose smallest mode is still evolving 
linearly. We used this small volume simulation for visual- 
ization and resolution studies. The 200 /i _1 Mpc simulation 
(mass resolution of 4.3 x 10 9 /i _1 Mg) is used for computing 
the quantitative results since a larger volume gives better 
statistics. 

The simulations were performed using the public ver- 
sion of the parallel Tree-PM code Gadget2 (Springel 2005) 



on a Linux cluster at the University of Groningen, Nether- 
lands. The initial conditions for both simulations were gen- 
erated at the z — 50 redshift using the transfer function 
given by Bardeen et al. (1986). 

5.1 Density and velocity divergence fields 

The output of the N-body simulation consists of a discrete 
set of particles. This needs to be interpolated to a continuum 
volume-filling density and velocity divergence fields that will 
be used as input for the Cosmic Web environment detection 
algorithm. It is crucial for the environment detection pro- 
cedure, especially for anisotropic features such as filaments 
and walls, that the interpolation method used to obtain the 
continuous fields retains all the scale and geometry informa- 
tion of the discrete galaxy or particle distribution. 

For these reasons we use the Delaunay Tessellation Field 
Estimator (DTFE), introduced by Schaap & van de Wey- 
gaert (2000) (for additional details see van de Weygaert & 
Schaap 2009; Cautun & van de Weygaert 2011), to recon- 
struct the underlying density and velocity divergence fields. 
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Figure 9. A 3D volume rendering of the density l + <5 (left panel) 
and velocity divergence 6 (right panel). Note that we only show 
the negative values of the velocity divergence. The picture rep- 
resents a 100 X 100 X 10( /i -1 Mpc) 3 volume in an N-body sim- 
ulation. We used the same volume to illustrate the Cosmic Web 
environments in Figures 12-16. 

For the environment detection algorithm, the DTFE method 
has the following important advantages: 

• Preserves the multi-scale character of the discrete dis- 
tribution. 

• Preserves the local geometry of the discrete distribu- 
tion. 

• Does not depend on user defined parameters or choices. 

The continuous DTFE density and velocity diver- 
gence fields are sampled on a 256 3 and 512 3 grid for the 
100 /i _1 Mpc and 200 /i _1 Mpc simulations respectively, such 
that there is a 0.4 /t -1 Mpc grid spacing in both cases. Figure 
8 shows a thin slice of the grid sampled density and velocity 
divergence fields. Note the level of detail in the structures, 
even inside voids, and the one-to-one correspondence be- 
tween the DTFE density and velocity divergence features. 
Figure 9 gives a 3D rendering of the density and velocity 
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Figure 10. Comparison between the mass of cluster-like 
objects detected using the NEXUS_den, NEXUS_tidal and 
NEXUS+ methods. The gray vertical line delineates the objects 
with mass larger than 5 ■ 10 13 ft _1 M0 while the diagonal black 
line shows a one-to-one relationship. 

divergence fields in a larger volume - the same volume that 
later on will be used to visualize the Cosmic Web environ- 
ments. 

5.2 Halo and subhalo catalogues 

We use the AMIGAs Halo Finder (AHF) by Knollmann & 
Knebe (2009) to identify the dark matter haloes. The AHF 
halo finder is the successor of the MHF halo finder by Gill 
ct al. (2004) . AHF uses adaptive mesh refinement to identify 
the density peaks which it classifies as the halo and subhalo 
centres. Afterwards it grows the objects around their centres 
until the spherically averaged density contrast reaches the 
virial density 5 . The last step consist in removing the grav- 
itationally unbound particles. The AHF halo and subhalo 
catalogues are complete up to haloes with 50 or more par- 
ticles (for a complete description see Knollmann & Knebe 
2009). 



6 THE COSMIC WEB ENVIRONMENTS 
6.1 Clusters 

The point-like objects detected by the NEXUS and 
NEXUS+ methods correspond to large overdensities in the 
density field. They range from very massive to very small 
mass objects. We will see later in this section that these ob- 
jects correspond to dark matter haloes. On the other hand, 
the cosmic clusters are the largest and most recent to form 
virialized objects (Voit 2005). To be able to identify our 
point-like objects with actual clusters, we need to limit our 
detections to only the most massive objects. In this study 
we consider as clusters the objects with mass larger than 
5 • 10 13 /i _1 Mg. This is a compromise between studying the 
most massive objects and having a large sample of such ob- 
jects in our simulation. 

5 The virial density is automatically computed by AHF and de- 
pends on both cosmological parameters and redshift. 
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Figure 11. Comparison of cluster-like objects with their corre- 
sponding AHF haloes. The top panel shows the mass compari- 
son while the middle panel shows the volume comparison for the 
NEXUS+ and NEXUS_den objects. The gray vertical line delin- 
eates the objects with mass larger than 5 • 10 13 h~ 1 M© while the 
diagonal black line shows a one-to-one relationship. The lower 
panel gives the mismatch in center position and mass difference 
between the AHF haloes and NEXUS_den clusters (only for ob- 
jects with mass larger than 5 ■ 10 13 /i -1 Mg). 



The cluster detection method performs very well for 
the NEXUS_den, NEXUS_tidal and NEXUS+ methods. 
In contrast, it does not give reliable detections for the 
NEXUS_veldiv and NEXUS _velshear methods. We suspect 
this deficiency is due to high vorticity in the cluster regions, 
vorticity that is not captured in the velocity divergence or 



shear fields and therefore is not taken into account in our 
methods. 

The NEXUS.den, NEXUS.tidal and NEXUS + methods 
detect the same clusters, with similar mass and volume as- 
sociated to each object. This can clearly be seen in the mass 
comparison plot shown in Figure 10. The NEXUS.den and 
NEXUS_tidal methods give very similar results, while we 
find a larger scatter when comparing these results with 
NEXUS+. 

There is a very close connection between the cluster-like 
objects of our methods and the most massive dark matter 
haloes. This is shown in Figure 11 where we use AHF to iden- 
tify the dark matter haloes. The top panel shows that there 
is a very good correlation between the AHF halo mass and 
its corresponding object identified using our methods. We 
show that this relation holds down to masses of 10 13 /i _1 M . 
This lower limit is a limitation of the grid spacing size 
and not of the method. We see that NEXUS.den, and also 
NEXUS_tidal (not shown), give a much better correlation 
of the cluster mass with the corresponding AHF halo mass, 
while in the case of NEXUS+ there is a larger scatter. A 
similar comparison is done between the cluster volumes and 
the AHF halo volume - see center frame of Figure 1 1 . Again 
there is a good match between the methods, but with a much 
wider scatter than in the case of the mass comparison. This 
is understandable since the volume is much more sensitive 
to the lower density regions around clusters, while the mass 
is dominated by the inner high density regions. In this case 
the scatter is especially large in the NEXUS+ results, than 
for the NEXUS_den and NEXUS_tidal methods. 

Finally we compare the NEXUS.den results and AHF 
haloes only for the most massive objects. In the bottom 
panel of Figure 11 it can be seen that the clusters agree 
within 10% to their corresponding AHF halo mass and that 
their centres are at most 0.4 grid spacing distance from 
the AHF halo center. The few outliers are objects that are 
merged in one of the methods and detected as distinct ob- 
jects in the second one. 

We see that NEXUS+ is less reliable in the detec- 
tion of Cosmic Web clusters than the NEXUS.den or 
NEXUS_tidal methods. The very localized nature of the Log- 
Gaussian filter means that the contribution of the highly 
dense center does not have a large effect on the periph- 
ery of the cluster. Because of this the outer boundaries of 
the NEXUS + clusters are more dependent on the substruc- 
ture at the periphery. One can overcome this "weakness" 
in the NEXUS+ method by identifying the clusters using 
NEXUS_den or NEXUS.tidal. 

6.2 Filaments 

After clusters, filaments are the most noticeable feature of 
the Cosmic Web. These structures are shown in Figure 12, 
as identified by the six methods proposed in this work. To 
obtain these results we have restricted our discussion to the 
most significant filaments, which we define as any continuous 
region with a volume larger than 10 (/i _1 Mpc) 3 . By doing 
so we discard small regions, typically around isolated haloes 
which even though show a local filament signature, are not 
embedded in the larger network. 

Comparing the filamentary maps with the density field 
rendering of the same volume from Figure 9, we see that 
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a) NEXUS_den d) NEXUS veldiv 




Figure 12. A 3D rendering of the filaments in a 100 X 100 X 10 (h 1 Mpc) 3 volume of the N-body simulation. The faint background 
shows the density field. The filaments were obtained using: a) NEXUS_den, b) NEXUS.tidal, c) NEXUS_denlog, d) NEXUS_veldiv, e) 
NEXUS_velshear and f) NEXUS+. This is the same region of the simulation volume as the density field shown in Figure 9. 
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all the methods succeed in identifying the strongest fila- 
ments, but there are differences when it comes to the smaller, 
less pronounced filaments. We immediately observe that the 
Cosmic Web filaments form an interconnected network, with 
the most massive filaments acting as the backbone. These 
very pronounced filaments branch into thinner ones that 
slowly disappear into lower density regions. The backbone 
filaments are clearly visible in the central panel while the 
branching into fainter structures is most pronounced in the 
lower-right panel of Figure 12. A visual inspection of clusters 
(not shown here) shows that these reside at the intersection 
of the most prominent filaments, which serve as the high- 
ways along which matter is transported to the clusters. 

A visual comparison of the filamentary maps obtained 
using the six methods leads to the following conclusions: 

• The filamentary structures in the NEXUS_den and 
NEXUSjveldiv are more clumpy and have less large 
scale cohesion when compared to the NEXUS_tidal and 
NEXUS_velshear results. This is in line with our expecta- 
tions, since the large scale modes bring a larger contribution 
for the latter methods (see Figure 6). 

• The methods using the tidal and velocity shear fields 
are biased towards the most significant structures and miss 
most of the filaments present in the less dense regions. 
This less pronounced filaments, marginally seen by the 
NEXUS_den and NEXUS_veldiv methods, are very well re- 
produced in the NEXUS_denlog and NEXUS+ methods. 
This supports the view that the Cosmic Web has structures 
present over a large range in density values. 

• While the most pronounced filaments are detected by 
all the methods, their thickness varies between the different 
methods. The NEXUS_denlog and NEXUS+ filaments have 
typical diameters around 2 h,~ 1 Mpc, while for the rest the 
typical diameter is around 4 /i _1 Mpc. The thinner filaments 
mostly constitute the inner regions of the much thicker fila- 
ments detected using the other methods. 

• While NEXUS_denlog finds the same structures as the 
other methods in the high density regions, it finds a much 
richer filamentary network in the underdense regions. This 
is because any small changes in the density field in these 
regions can lead to a large contrast in the density loga- 
rithm. The Poissonian sampling noise is especially important 
for density determination in the underdense regions due to 
sparse sampling. This makes the NEXUS_denlog method es- 
pecially sensitive to Poissonian noise in the void-like regions. 

To summarize, the variation in the six methods mani- 
fests itself as mostly differences in the detection of smaller 
filaments and thickness differences in the very prominent 
structures. Since the large scale modes contribute much 
more to NEXUS_tidal and NEXUS.velshear, these meth- 
ods identify only the largest filaments that correspond to 
the peaks of the large scale modes. On the other hand, 
the NEXUS_denlog and NEXUS+ methods are much more 
sensitive to the less pronounced structures, finding an im- 
portant filamentary network also in the underdense regions. 
While NEXUS_den and NEXUS.veldiv are in between the 
two classes of results, they are much closer in character to 
the NEXUS -tidal and NEXUS_velshear methods. 




Figure 13. The dark matter particles left after taking out the 
particles located in cluster and filament environments. This is a 
projection of the volume of the N-body simulation as in Figures 
9, 12 and 14. 

6.3 Walls 

The wall environments detected using the six methods are 
presented in Figure 14 where, on top of the walls in orange, 
we superimposed the filaments in light blue. In the case of 
filament identification one can use the density and veloc- 
ity divergence maps to judge the success of the detection 
method, but this is much more difficult for walls. This is 
a consequence of both the smaller contrast and the planar 
nature of these structures. The presence of sheet-like struc- 
tures in the distribution of matter on cosmological scales 
can be easily inferred from Figure 13, where we show the 
dark matter particles after removing the particles located in 
clusters and filaments. The most striking structures are the 
line-like arrangements of particles visible especially in the 
upper part of the figure. These are sheets that are perpen- 
dicular on the projection plane and not filaments missed by 
our detection algorithm. There are additional walls along 
the projection plane (e.g. center upper part, between the 
three line-like structures) but these are less easily detected 
visually. 

By comparing the particle distribution from Figure 13 
to the NEXUS walls in Figure 14 it is clear that the algo- 
rithm is very successful in identifying both the prominent 
as well as the tenuous Cosmic Web walls. The resulting ob- 
jects form large continuous planar structures which delineate 
the different voids. The aspect of the sheets resembles very 
much that of walls of biological cells that have been stacked 
on top of each other, with the cells having a wide range of 
sizes. Even though the Cosmic Web walls are continuous for 
large portions of their surface, they are still punctured by a 
large number of holes that allow for void percolation. While 
these holes are dependent on the wall identification method, 
they are a sign of the diffuse and sparse nature of the cosmic 
walls. 

As in the case of the filamentary network, the prominent 
walls are detected by all the methods, but there are differ- 
ences when it comes to the sheets in the more underdense 
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a) NEXUS den 



d) NEXUS_veldiv 




Figure 14. A 3D rendering of the walls (orange) in a 100 X 100 X 10 (/i _1 Mpc) 3 volume of the N-body simulation. The light blue depicts 
the filaments while the faint background shows the density field. The walls were obtained using: a) NEXUS-den, b) NEXUS-tidal, c) 
NEXUS_denlog, d) NEXUS_veldiv, e) NEXUS_velshear and f) NEXUS+. This is the same region of the simulation volume as the density 
field shown in Figure 9 and the filaments in Figure 12. 
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Figure 15. The dark matter particles in the different environments of the Cosmic Web. The panels gives: a) all, b) filament-only, c) 
wall-only and d) void-only particles. The environments where identified using the NEXUS+ method. 



regions. Some of the most important features and differences 
between the six results can be summarized as: 



• The NEXUS_den walls have a very clumpy appear- 
ance and this is also true to a lesser extent for the 
NEXUSTidal results. This is in contrast with the other 
methods where the walls have a much more planar and sheet- 
like look. The clumpy appearance is due to the composition 
of walls, which are tenuous structures with sporadic haloes 
from place to place. The concentration of mass in the haloes 
compared to their neighbourhood regions gives the clumpy 
structure of the NEXUS.den and NEXUSTidal walls. 

• The NEXUS_veldiv and NEXUS _velshear walls have a 
smooth planar appearance since the velocity field, compared 
to the density, is less affected by the sparse Poissonian sam- 
pling of the mildly underdense regions. These are the re- 
gions that make most of the volume in walls. We find the 



same smooth planar look also for the NEXUS_denlog and 
NEXUS+ walls. 

• We find that for all the methods the larger gaps in the 
walls are present at the same locations, indicating that, at 
least occasionally, there is no clear boundary between adja- 
cent voids. 

• The NEXUS _velshear method seems to be the most con- 
servative tracer of walls. It detects only the most prominent 
wall regions. 

A visual inspection favours the NEXUS+, 
NEXUS_denlog and NEXUS_veldiv as best tracers of 
the cosmic sheets with the remaining methods either giving 
clumpy detections or missing some of the more tenuous 
structures. 

By comparing the filaments and walls in Figure 13 we 
find that most if not all the filament volume elements are 
embedded in walls. Though not shown, we find that also 
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Figure 16. The dark matter haloes as a function of the environment they reside in. The panels gives: a) all, b) filament-only, c) wall-only 
and d) void-only haloes. The environments where identified using the NEXUS+ method. The colour and size of the points is proportional 
to the mass of the halo they represent. The labels in the legend correspond to log 10 (M/ /i —1 Mq). Note that most of the smaller mass 
haloes around massive ones are not visible due to the larger size with which we show the massive haloes. 



the clusters are fully embedded in filaments. This suggests 
that the strict classification where each volume element is 
assigned to a single environment can be extended by realiz- 
ing that clusters are embedded in filaments and that in turn 
filaments are embedded in walls. 



6.4 Cosmic Web dark matter particle and halo 
population 

An interesting and important issue for the study of struc- 
ture formation is our understanding on how far the cosmic 
structure in the dark matter distribution is reflected in the 
halo distribution. In Figure 15 we show the particles in our 
numerical simulation splitted according to the environment 
in which they reside. The filament environment correctly 
traces the largest linear particle concentrations. This can 
clearly be seen in the prominent filaments. While the thin- 



ner structures are less populated with dark matter parti- 
cles, they are identified as filaments due to their higher lo- 
cal contrast. NEXUS+ finds that on average filaments have 
an overdensity 1 + 6 around 11, as can be seen by compar- 
ing the mass and volume filling fractions given in Table 3. 
On the other hand, walls are much more tenuous structures 
and this is visible in the particle distribution. When seen 
edge on, the sheets appear as prominent structures in the 
particle distribution. But when looked at face on, they are 
sparsely populated by particles and hence difficult to iden- 
tify. When comparing the particle distributions in walls and 
voids there seems to be little difference in the particle den- 
sities. This is just a projection effect, with walls having on 
average an overdensity around 1.4 versus an overdensity of 
0.2 for voids (see Table 3 for details). The success of the 
structure finding method can be easily assessed by observ- 
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Table 3. The mass and volume filling fractions for the environ- 
ments identified using NEXUS+ . The results presented here are 
for the larger 200 h~ 1 Mpc simulation. The mass and volume fill- 
ing fractions for the other five methods are given in appendix 
A. 



Environment 


Mass fraction (%) 


Volume fraction (%) 


clusters 


8.0 


0.027 


filaments 


51.3 


4.35 


walls 


24.0 


16.8 


voids 


16.7 


78.8 



10" F 




M(M 8 h-') 



Figure 17. The cumulative halo mass function segmented ac- 
cording to the components of the Cosmic Web as identified by 
NEXUS+. The upper panel gives the cumulative mass function 
normalized according to the volume of the whole simulation box. 
The lower panel gives the cumulative mass function normalized 
according to the volume of each environment (see Table 3). 



ing that the dark matter particles in void regions do not 
have any significant structure present in their distribution. 

A more interesting picture is found when looking at the 
variation of the dark matter halo populations with the Cos- 
mic Web environment. This is shown in Figure 16 where 
the haloes are coloured and scaled according to their mass. 
Figure 16, in conjunction with Figure 17 which gives the 
cumulative halo mass function split according to environ- 
ment, offers a very suggestive picture. We find that all of 
the massive haloes with M ^ 5 x 10 13 /i _1 M0 are located 



in cluster environments . Moreover, clusters are the most 
crowded regions when it comes to haloes of all masses, with 
halo overdensities about 10 times larger than in filaments 
and about 100 times larger than in the full simulation box 
(see lower panel in Figure 17). 

When looking at filaments we see that these environ- 
ments are also crowded when it comes to haloes. Most of 
the 10 12 /i _1 Mq and higher mass haloes are located in fila- 
ments, while the lower mass objects are a factor of 10 more 
common in filaments than on average in the universe. From 
Figure 16 we see that even the more tenuous filaments have a 
large number of haloes which shows the power of our method 
to correctly identify the filamentary environments. 

The walls are dominated by low mass objects, with the 
sheets containing only a significant share of the 10 12 /i _1 Mg 
and lower mass haloes. For a few times 10 11 /i -1 M0 and 
lower mass haloes, the walls have a similar halo density as 
the average universe. Compared to voids, the sheets clearly 
have a much higher halo density and are populated by more 
massive haloes. In contrast, the voids are very sparsely pop- 
ulated, with extremely few 10 11 /i _1 Mq and higher mass 
haloes. As in the case of the particle distribution, there does 
not seem to be any significant structures present in the void 
halo distribution. 



6.5 Single scale versus multiscale analysis 

One of the frequent questions that surface when dealing with 
cosmological structure identification is the optimal value of 
the smoothing scale (see Hahn et al. 2007a; Forero-Romero 
et al. 2009). Often this constitutes a major limitation of 
single scale approaches. Moreover, until now, the choice of 
one smoothing length versus another was mostly heuristic. 
While our multiscale approach does not suffer from these 
problems, it is certainly a very interesting question to find 
which scales are the most important and how they compare 
with results from other works. 

We exemplify the effects of a single scale versus a multi- 
scale approach in Figure 18. To obtain the single scale results 
we restrict the NEXUS+ method to a single filter radius. We 
then show this result for several values of the filter radius. 
For small smoothing scales, the method detects most of the 
filaments. While it reproduces correctly the small filaments, 
it greatly underestimates the thickness of the more promi- 
nent objects. As the filter radius is increased, most of the 
thinner filaments are missed while the larger ones are de- 
tected as being thicker and thicker. In conclusion, applying 
a filter of a given radius makes the environment detection 
method sensitive only to objects similar or larger than that 
smoothing scale. Moreover, at a given filter radius, all fila- 
ments have very similar diameters which are given by the 
value of the smoothing radius. 

By comparing the single-scale results versus the full 
NEXUS+ results (see Figure 18 panel d) we conclude that 
the 2 /i _1 Mpc smoothing scale gives the closest match to the 

6 Remember that the 5 X 10 13 /i -1 M0 mass threshold was in- 
troduced as a lower mass cut-off in section 6.1. While this lower 
mass threshold depends on the cluster definition one chooses, the 
method is very successful in identifying as clusters all the haloes 
above the cut-off mass. 
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a> 1 Mpc/h 




c) 4 Mpc/h 



d) NEXUS+ 




Figure 18. A comparison of the single-scale versus the multi- 
scale approach. Panels a) to c) show the single scale filaments 
obtained by restricting the NEXUS+ method to the single filter 
radius of 1, 2 and 4 h~ 1 tApc respectively. Panel d) shows the full 
NEXUS+ results. 



multiscale filaments. However, there are a lot of important 
differences, with many thin filaments missing in the single 
scale picture. Another striking difference is the smaller di- 
ameters for the larger filaments in the single scale versus 
the full NEXUS+ results. These results strongly show the 
need of a multiscale approach to be able to fully trace all 
the features of the Cosmic Web. 

When restricted to a single scale, the NEXUS_den and 
NEXUS_veldiv methods show an even stronger difference be- 
tween results at different filter radii. This trend is so strong 
that it is very difficult to find a single smoothing scale that 
matches even remotely the multiscale results. Even the sin- 
gle scale NEXUS_tidal and NEXUS_velshear methods show 
significant differences between scales, though less prominent 
than the results in Figure 18. The weaker dependence on 
scale is due to the steep drop at small wavelengths of the 
gravitational and velocity potentials spectra (see Figure 6) 
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Figure 19. The NEXUS+ filaments coloured according to the 
smoothing scale that gives the largest filamentary signature for 
that voxel. The corresponding colour bar indicates the filter scale 
in units of h~ Mpc. The lower panels b) to d) show a smaller 
region from the larger volume presented in panel a). These panels 
show the filaments detected using a maximum of: b) 1 /i _1 Mpc, 
c) 2 h~ iJVIpc and c) 4 h~ iJVIpc smoothing scales. 



which means that large scale modes contribute more than 
for the other methods. For the single scale NEXUS_tidal and 
NEXUSjvelshear methods we find that 1 /i _1 Mpc smooth- 
ing radius results are the closest match to the full multiscale 
results. This is in contrast with Hahn et al. (2007a), who ar- 
gue for a 2.1 /i _1 Mpc filter radius. 

Another very interesting question is finding the smooth- 
ing scale at which a region has the largest environmental sig- 
nature. This is illustrated in Figure 19. We show in different 
colours the filter radius which gives the strongest filamen- 
tary characteristic for a given region. The small filter radii 
give the strongest signature in the thinner filaments while 
the larger smoothing scales give a stronger signal for the ma- 
jor filaments. A closer inspection of the prominent filaments 
shows that their central axes give the strongest filamentary 
response for small filter radii. Increasing the size of the filter 
adds larger and larger filamentary regions around this inner 
central axis. This is visible in panels b) to d) of Figure 19. 
This is due to the inner structure of the filaments, with their 
central axis having a larger density than their periphery. 



7 COMPARISON TO OTHER STRUCTURE 
FINDING ALGORITHMS 

In this section we discuss the abilities and virtues of the 
NEXUS and NEXUS + algorithms with respect to those 
of other structure finding algorithms. It will underline the 
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Figure 20. The filamentary network as identified by the al- 
gorithms presented in: a) Halm et al. (2007a) and b) Zhang 
et al. (2009) (see text for details). The results are obtained for a 
2 /i _1 Mpc filter width. It shows the same volume of the N-body 
simulation as in Figures 9 and 12. 



advantages and disadvantages of the instruments described 
and introduced in this paper. In our discussion we limit our- 
selves to Hessian and topological based methods. There are 
a few additional methods that find filaments using directly 
the particle/galaxy distribution (Stoica et al. 2005, 2007, 
2010; Chazal et al. 2009). While these approaches have the 
advantage that one does not need to compute the density 
field, they depend on many free parameters that make their 
use cumbersome. 

A first group of Cosmic Web identification methods con- 
sists of algorithms that use the Hessian of the density or 
gravitational potential. Both the NEXUS and NEXUS+ al- 
gorithms are part of this class. 

Hahn et al. (2007a) (referred to as HPCD) proposed the 
use of the tidal field eigenvalues for environment classifica- 



tion . The method uses the criterion that: all positive eigen- 
values identify clusters, one negative eigenvalue corresponds 
to filaments while two negative eigenvalues trace sheets. The 
result of this method applied to our simulation is shown in 
Figure 20. It is immediately clear that while the inner region 
of the HPCD results correspond to NEXUS_tidal (compare 
to panel b in Figure 12), the HPCD filaments extend to much 
larger diameters, encompassing substantial parts of wall and 
void regions. This leads to a large cross-contamination of the 
Cosmic Web components. A similar method to HPCD , but 
using the density Hessian instead of the tidal field, was used 
by Zhang et al. (2009). From Figure 20, panel b, is obvi- 
ous that this last method is not very successful. It leads 
to even more misclassified regions than the HPCD method. 
There are two main shortcomings of the two methods: use 
of a single scale approach and the absence of a threshold to 
distinguish between significant and spurious detections. 

The absence of a detection threshold in the 
HPCD method was pointed out in Forero-Romero et al. 
(2009). They suggested that an eigenvalue threshold in the 
range 0.2 — 0.4 gives results in agreement with the visual im- 
pression of the Cosmic Web. Using NEXUS_tidal restricted 
to the 2 /i -1 Mpc filter radius we do obtain that the filament 
and wall threshold is 0.42 and 0.2 respectively. While we 
do confirm the results of Forero-Romero et al. (2009), we 
stress that NEXUS_tidal has the advantage of a multiscale 
approach and does not employ user dependent arguments 
for specifying the detection threshold. 

It is the multiscale character of NEXUS and 
NEXUS+ which are crucial for their successful analysis of 
emerging structures and patterns in the hierarchically evolv- 
ing cosmic mass distribution. The first version of the mul- 
tiscale formalism that we have developed into NEXUS and 
NEXUS+ is the Multiscale Morphology Filter (MMF), de- 
scribed and introduced in Aragon-Calvo et al. (2007b). 
NEXUS is an extension of MMF to a more physical and 
versatile algorithm. While MMF is restricted to the density 
field, the NEXUS method incorporates, among others, the 
tidal and velocity fields. Both MMF and NEXUS_den find 
the same filamentary structures, while the second method 
gives a much better identification of walls (for details see 
Appendix A). More importantly, NEXUS+ seems optimal 
at capturing the structural intricacies of the Cosmic Web, 
and thus represents an important advancement within the 
context of the multiscale scale-space formalism that we have 
been developing. 

An additional important class of structure identification 
procedures are based on topological considerations, in prin- 
cipal following an analysis of the Morse-Smale complex of 
the density field (Sousbie et al. 2008; Sousbie 2011; Aragon- 
Calvo et al. 2010). Given the fact that in general, cosmic 
density fields behave like a proper Morse function, the as- 
sumption is that filaments and walls in the mass distribution 
can be identified with the manifolds in the density field con- 
necting maxima via saddle points of the field Sousbie (2011). 
Filaments are identified with the line connecting two max- 
ima via a saddle point, and walls with the sheet separat- 
ing the regions around two minima and centered around a 



7 This method can be easily implemented as a special case of our 
NEXUS_tidal algorithm. 
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saddle point. Aragon-Calvo et al. (2010) developed a simi- 
lar strategy by first delineating the watershed basins around 
the mimima in the field, following the WVF procedure intro- 
duced by Platen et al. (2007) in which the watershed basins 
are identified with the voids in the cosmic mass distribution. 
The topological character of the boundary is determined lo- 
cally via the number of touching watershed basins. Walls 
are the 2-d manifolds separating two watershed basins, the 
filaments are their boundaries and identified as the locations 
where three basins touch each other. 

NEXUS and NEXUS+ do find that the largest filaments 
are between massive density peaks, which is in agreement 
with the implicit assumptions of the topological methods. 
But on top of the prominent structures, we also find, es- 
pecially in the NEXUS+ results, an important network of 
thinner objects which branches into voids. These tenuous 
objects contradict the hypothesis that filaments are always 
located between density maxima. We also find that walls are 
not fully continuous sheets and that they sometimes stop as 
they branch into lower density regions. Thus, the voids fully 
percolate, with large regions without a clear boundary be- 
tween adjacent voids. Also, we should note that while the 
topological based methods can detect the central axis of fila- 
ments and inner plane of walls, they cannot assign a natural 
thickness to the structures. As our results show, both the 
filament and wall environments have a wide range of sizes. 
This makes it rather challenging for the topological methods 
to detect and outline in a natural fashion the regions belong- 
ing to a filament around the identified central filament axis. 



8 CONCLUSIONS AND PROSPECTS 

This work presents the NEXUS and NEXUS+ methods, 
which are multiscale and automatic algorithms used for the 
segmentation of the Cosmic Web into its distinct compo- 
nents: clusters, filaments, walls and voids. We have shown 
that the environments identified with the two methods cor- 
respond very well to the structures visible in the density 
and velocity divergence fields as well as in the dark matter 
particle and halo distributions. The success of the method 
lies in two important ingredients: the use of a scale free ap- 
proach that makes sure that the algorithm detects structures 
of all sizes, and the use of a physically motivated threshold 
to distinguish valid environments from spurious detections. 
Another strength of the two algorithms is that they do not 
depend on user set parameters and therefore can be easily 
applied in a consistent way to multiple data sets. 

We have extended the NEXUS method to detect the 
Cosmic Web as traced by the density (NEXUS_den), tidal 
field (NEXUS_tidal), velocity divergence (NEXUS_veldiv) 
and velocity shear fields (NEXUS_velshear). We find that 
NEXUS_den and NEXUS_tidal are very efficient in identify- 
ing the cluster regions when compared with the most mas- 
sive haloes present in the simulation. The methods perform 
similarly in the detection of filament and wall regions, with 
all the prominent structures detected by all methods. The 
only differences arise in the identification of the more ten- 
uous structures. We find that NEXUS+ performs better in 
tracing the weaker filaments and walls. 

From all the methods presented in this paper, we find 
that NEXUS+ is the most successful one in tracing the Cos- 



mic Web. Its main advantage comes from the use of the Log- 
Gaussian filter which is designed to better deal with the or- 
ders of magnitude difference in the density field between low 
and high density regions. The filamentary and wall like envi- 
ronments detected with NEXUS+ contain complex networks 
of prominent structures that branch out into more tenuous 
ones until they finally disappear out into underdense regions. 

We showed in section 7 that NEXUS and 
NEXUS+ have several advantages compared to other 
Hessian based methods: the use of a multiscale approach 
and a physically motivated threshold for identifying the 
significant environments. Compared to topological methods, 
our tools are able to detect the filamentary /wall regions 
and not only their central axis/plane. Moreover, we do 
not make the assumptions that filaments extend between 
density maxima and walls separate density minima basins. 
Our results show that these assumptions do not always 
hold. 

Equipped with the NEXUS and NEXUS+ methods, we 
plan to address a range of cosmological issues. The abil- 
ity to identify filaments and walls over a range of scales, 
in both numerical and observational datasets, allows us to 
study not only environmental factors affecting the formation 
and evolution of dark matter haloes and galaxies, but also 
the hierarchical buildup of the Cosmic Web itself. Our first 
priority lies with a systematic study of environmental factors 
affecting the evolution of galaxies. The fact that we can iden- 
tify galaxies and haloes within finely outlined filaments and 
sheets will allow us to determine which physical character- 
istics are most sensitive to the environment, which galaxies 
and haloes are most sensitive to large scale influences, and 
to study the causes that give rise to this dependency. 

A major point of our interest is performing a system- 
atic comparison between the structures traced by the den- 
sity, velocity and tidal fields in order to understand which 
of these physical influences are most decisive in determining 
the global outline of the Cosmic Web. In addition to our 
current focus on the large scale structure of the dark matter 
distribution, we will also direct our study to the structure 
of the gaseous Cosmic Web. Comparison of the IGM with 
the dark matter structures in numerical simulations will be 
instrumental in understanding how the Cosmic Web can be 
traced both in the galaxy and cosmic gas distribution. This 
will be essential for relating the distribution of HI in the lo- 
cal Universe to the overall large scale structure found in the 
galaxy distribution (Popping & Braun 2011), and will help 
understand recent findings such as a small-scale HI filament 
in a void (Beygu et al., in prep.). 

Finally, the NEXUS and NEXUS+ procedures are per- 
fectly suited for a systematic appraisal of the structures 
found in maps produced by galaxy redshift surveys such as 
SDSS and 2MRS. 
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APPENDIX A: OPTIMAL FILAMENT AND 
WALL DETECTION 

This section deals with identifying the signature threshold 
used for the detection of the Cosmic Web filaments and 
walls. When applying the environment detection algorithm 
every region of space is assigned an environment signature 
5 (for definition see Eq. (7)). A large signature corresponds 
to very prominent structures while a zero or small one cor- 
responds to null detections. Since many regions of space will 
have a signature value between the two extremes, we need 
to identify a signature threshold that differentiates between 
valid structures and spurious detections. All regions with 
signatures larger than the threshold will correspond to valid 
environments. 

The simplest way to define the Cosmic Web is using the 
tidal field, since this is what drives the anisotropic collapse. 
Using the eigenvalues of the tidal tensor one defines fila- 
ments and walls as regions with one and two negative eigen- 
values (see Hahn et al. 2007a). The major problem with this 
approach is that it gives unrealistic looking environments, 
with only a small fraction of the volume (~ 10%) occupied 
by voids - this is in stark contrast with the observational 
data. 




Figure Al. A 3D volume rendering of the NEXUS_dcn filaments 
(blue) and walls (orange) when the detection threshold is taken 
as the percolation point. 
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Figure A2. The dependence of the mass fraction in the compo- 
nents of the Cosmic Web as a function of environment signature. 
The two panels shows the dependence for filaments and respec- 
tively walls obtained using the NEXUS_den method. The contin- 
uous red curve gives the mass fraction M in filaments/ walls while 
the dashed blue curve gives AM 2 (see text for definition). The 
black curve gives the volume fraction of the largest filament/wall. 
The sharp transition of this curve from to 1 marks the percola- 
tion threshold. The vertical gray lines mark the percolation point 
(filled line) and the peak of AAf 2 (dashed line). 
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Figure A3. Comparison of the mass fraction (upper panel) and 
volume fraction (lower panel) in filaments detected using the 6 
methods described in this paper. The continuous black line de- 
scribes the filaments detected using the AM 2 threshold method 
(described in this section) while the dashed red line shows the fil- 
aments detected using the percolation threshold (Aragon-Calvo 
ct al. 2007b). The horizontal dotted line give the average mass 
and volume fraction for the AM 2 results. 



AM = 



dM 2 



dlogS 



(Al) 



This is represented by the dashed blue curve in Figure A2. 
We found that the peak position in AM 2 gives a robust 
and natural way of identifying the most significant filaments 
and walls. The environments detected using the AM 2 peak 
threshold are shown in Figures 12 and 14. This new thresh- 
old captures very well the filamentary and wall features seen 
in the density and velocity divergence fields from Figure 9. 
Moreover the peak of AM 2 for the NEXUS_den filaments 
is very close to the percolation threshold and hence repro- 
duces the success of using percolation as a good tracer of 
the filamentary network 8 . 

There are two major improvements when using the 
AM 2 peak versus the percolation point as the detection 
threshold for the most significant environments. The walls 
detected via the AM 2 method have a bigger large scale co- 
hesion than the results of the percolation method. This can 
clearly be seen by comparing the NEXUS_den walls from 
Figure 14 and the percolation walls in Figure Al. The second 
enhancement comes from a more robust detection threshold. 
This can be seen in Figure A3 where we compare the mass 
and volume fraction in filaments detected using the 6 meth- 
ods introduced in this paper. Notice the large variation in 
mass fraction between the NEXUS_den and NEXUS+ re- 
sults for the percolation method - with the former having 
almost twice as much mass than the latter one. In the case of 
the AM 2 threshold all methods give a similar mass fraction, 
with a much smaller scatter around the mean. The same be- 
haviour can be seen in the volume fraction plot, where the 
AM 2 threshold gives more consistent values. The only ex- 
ception is the NEXUS + volume fraction whose lower value 
is due to the method itself and not the detection threshold 
used (see section 6.2 for more details). 



A second approach is to detect only the most signifi- 
cant filaments and walls. This is the procedure we follow in 
this paper. Aragon-Calvo et al. (2007b) have argued that 
the percolation threshold of filaments/ walls offers a natural 
way of identifying the prominent structures. We find that 
indeed the percolation threshold gives a good identification 
of the filamentary network, but fails in detecting the sheets. 
The walls obtained via this method are made of many small 
patches that do not show the large scale cohesion expected 
for void boundaries (see Figure Al). Another downside is 
the dependence of the percolation threshold on the grid res- 
olution used to analyse the data. But more importantly for 
this work, many large differences between the results of the 
6 methods described here can be attributed to different per- 
colation properties of the filamentary /wall networks in each 
method. We analyse this in more details later on. 

Most of the mass in filaments and walls is given by re- 
gions with a narrow range in environment signatures. This 
is illustrated by the red curve of Figure A2 which gives the 
mass fraction M in filament/wall regions as a function of 
the signature threshold S. The rapid increase in the fil- 
ament/wall mass can be appreciated when computing the 



8 The percolation threshold and the AM 2 peak are similar only 
for the NEXUS_den filaments, in general there are large offsets 
between the two values. This was expected since Aragon-Calvo 
et al. (2007b) have shown the success of the percolation threshold 
for filament detection only when using the density field as tracer 
of Cosmic Web environments. 



